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This  paper  introduces  an  optimization  method  for  improving  thermoelectric  heat  pump  performance  by 
operating  condition  management  of  the  thermoelectric  modules  (TEMs)  and  design  optimization  of  the 
heat  exchangers  linked  to  the  TEMs.  The  device  studied,  corresponding  to  an  original  configuration  of 
the  thermoelectric  heat  pump,  comprises  two  commercial  thermoelectric  modules  and  two  mini-channel 
heat  sinks  through  which  water  flows,  in  contact  with  both  sides  of  the  TEMs.  The  objective  function  is 
the  maximization  of  the  device’s  coefficient  of  performance  (COP),  including  the  electrical  and  mechan¬ 
ical  consumption  of  the  thermoelectric  modules  and  the  circulating  auxiliaries.  First,  the  optimization 
variables  are  the  number  and  the  diameter  of  mini-channels,  and  the  mass  flows  for  both  heat  sinks 
(hot  and  cold  sides).  The  results  show  that  similar  results  are  obtained  by  minimization  of  the  entropy 
generation  in  the  device.  Finally,  the  hot  thermal  power  demand  is  included  in  the  optimization  variables 
for  complete  optimization  of  the  device.  The  results  of  full  optimization  converge  with  those  obtained 
with  the  previous  partial  optimization. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  heat  pumps  are  devices  that  exploit  the  temper¬ 
ature  difference  created  by  the  Peltier  effect  between  the  junctions 
of  semi-conductors.  Thermoelectric  phenomena  [1-3]  were  ob¬ 
served  and  described  during  the  first  half  of  the  19th  century. 
However,  because  of  the  low  figure  of  merit  characterizing  the 
quality  of  the  thermoelectric  material  used  at  this  time,  the  ther¬ 
moelectric  effect  was  mainly  exploited  for  measuring  high  temper¬ 
atures.  In  1885,  Rayleigh  proposed  creating  a  thermoelectric 
generator  [4],  but  the  first  real  application  only  appeared  in  1947 
with  the  Telkes  thermoelectric  generator  [5].  During  the  period 
from  1957  to  1965,  many  thermoelectric  materials  were  discov¬ 
ered,  notably  bismuth  telluride,  which  shows  the  highest  figure 
of  merit  at  ambient  temperature  [6].  Today,  this  material  is  still 
the  most  widely  used  for  applications  at  ambient  temperature.  In 
a  thermoelectric  heat  pump,  the  two  important  characteristics 
are  the  quality  of  the  material  used  and  the  thermal  resistances 
of  the  device.  Many  studies  have  been  conducted  on  the  optimiza¬ 
tion  of  the  thermoelectric  device;  however,  most  of  them  deal  with 
only  one  part  of  the  device,  such  as  the  thermoelectric  elements 
[7-9],  the  contact  thermal  resistance  [7],  or  the  influence  of  the 
heat  exchangers  [10,11].  Crane  and  Jackson  [12]  also  worked  on 
the  optimization  of  cross-flow  heat  exchangers  coupled  to  a  ther¬ 
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moelectric  module  with  the  ratio  between  the  useful  thermal 
power  and  the  cost  of  the  device  as  the  objective  function.  Khire 
et  al.  [13]  worked  on  the  optimization  of  a  thermoelectric  heat 
pump  coupled  to  a  building.  The  optimization  simultaneously  min¬ 
imizes  two  design  objectives,  the  input  power  and  the  number  of 
thermoelectric  modules.  In  this  study,  a  method  based  on  coupling 
the  thermoelectric  phenomena  to  heat  transfer  and  pressure  drop 
phenomena  in  the  heat  exchangers  is  proposed  to  optimize  a  ther¬ 
moelectric  heat  pump  operating  with  water.  This  method  allows 
precise  determination  of  the  optimal  geometrical  dimensions  of 
the  heat  exchangers  and  the  optimal  operating  conditions  for  given 
fluid  temperatures.  The  objective  function  considered  is  the  de¬ 
vice’s  coefficient  of  performance  (COP)  including  the  mechanical 
consumption  of  the  circulating  auxiliaries.  Design,  mass  flows 
and  electrical  power  are  considered  as  optimization  variables.  Fi¬ 
nally,  the  results  obtained  are  compared  with  a  second  objective 
function:  the  minimization  of  entropy  generation.  The  potential 
application  concerns  heating  and  cooling  in  low-consumption 
buildings. 

2.  Configuration  of  the  device 

The  thermoelectric  heat  pump  unit  considered  in  this  study 
comprises  three  heat  exchangers  and  two  thermoelectric  modules. 
This  thermoelectric  heat  pump  configuration  results  from  a  previ¬ 
ous  study  on  an  optimized  method  of  thermoelectric  heat  pump 
management  designed  to  increase  the  global  performance  of  the 
device  [14].  This  study  showed  that  it  is  possible  to  maintain  an 
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Nomenclature 

A 

surface  area,  m2 

Subscripts 

Cp 

specific  heat,  J  Kg-1  K_1 

c 

cold 

D 

diameter,  m 

elec 

electrical 

E 

space  between  two  channels,  m 

gen 

generate 

i 

friction  factor 

h 

hot 

h 

heat  transfer  coefficient,  W  m-2  K_1 

opt 

optimal 

H 

height,  m 

m 

mechanical 

I 

electrical  intensity,  A 

device 

relative  to  the  complete  device 

I< 

thermal  conductance,  W  K_1 

L 

length,  m 

Superscripts 

m 

mass  flow,  Kg  s-1 

b 

bulk 

N 

number  of  channels 

cer 

ceramic 

Nu 

Nusselt  number 

etc 

contact 

P 

pressure 

ch 

channel 

Pr 

Prandtl  number 

cond 

conductivity 

R 

resistance,  KW'1 

conv 

convection 

Re 

Reynolds  number 

f 

fluid 

S 

entropy,  W  K_1 

in 

inlet 

T 

temperature,  K 

int 

interface 

V 

velocity,  m  s-1 

j 

junction 

out 

outlet 

Greek  symbols 

tot 

total 

OC 

Seebeck  coefficient,  V  K_1 

opt 

optimal 

2 

thermal  conductivity,  W  m_1  K_1 

TEM 

thermoelectric  module 

* 

power,  W 

P 

density,  Kg  m-3 

optimal  COP  of  the  device  for  given  temperature  conditions  using 
parallel  and  cascade  associations  for  any  thermal  power  needed. 
The  configuration  presented  in  Fig.  1  illustrates  a  parallel  operating 
mode.  Heat  transfer  is  mainly  ensured  by  convection  in  the  work¬ 
ing  fluid.  In  this  configuration,  the  same  mass  flow  for  the  two  cold 
heat  exchangers  and  the  same  electrical  intensity  supplied  for  both 
TEMs  are  assumed.  Finally,  the  hot  heat  exchanger  dissipates  the 
same  thermal  power  on  its  two  sides. 

Fig.  2  illustrates  the  same  device  presented  in  Fig.  1  in  cascade 
operating  mode.  In  this  configuration,  the  two  TEMs  are  associated 
in  cascade  and  the  thermal  flux  is  transferred  from  the  cold  heat 
exchanger  to  the  hot  heat  exchanger  via  thermal  conduction 
through  the  central  heat  exchanger.  The  working  fluid  does  not 
flow  through  the  hot  heat  exchanger.  This  configuration  allows 
increasing  the  maximal  COP  with  a  reduction  of  the  corresponding 
thermal  power.  Consequently,  the  heat  exchanger  studied  here  is 
original  because  two  operating  modes  are  considered:  conduction 
through  the  bulk  when  the  thermoelectric  modules  are  associated 
in  cascade  and  convection  when  the  thermoelectric  modules  are 


USEFUL  SIDE  ( e  g.  building) 

300  K 

Fig.  1.  Scheme  of  the  device  for  a  parallel  mode. 


TEM 

Heat  exchanger 
Thermal  power  direction 

Hot  fluid 
Cold  fluid 
Electrical  network 


Fig.  2.  Scheme  of  the  device  for  a  cascade  mode. 


associated  in  parallel.  This  original  heat  exchanger  is  covered  by 
patent  FR  11  53394. 

Heat  transfer  by  conduction  in  heat  exchangers  (cascade)  is  less 
restrictive  than  by  convection  (parallel);  consequently,  this  study 
focuses  on  the  parallel  mode.  Therefore  the  method  introduced 
here  aims  to  optimize  the  design  of  the  heat  exchangers  and  the 
operating  conditions  for  the  parallel  operating  mode. 

Only  the  heating  mode  is  considered  here.  The  hot  and  cold 
fluid  temperatures  are  respectively  set  to  300  K  and  285  K 
(Fig.  1 ),  corresponding  to  floor  heating  temperature  in  the  building 
and  to  a  constant  geothermal  source,  respectively.  These  operating 
temperatures  are  kept  constant  in  the  following  in  order  to  validate 
the  optimization  approach. 

3.  Analytical  model 

The  goal  of  this  study  is  to  optimize  the  performance  of  a 
thermoelectric  heat  pump  device.  The  optimization  method  must 
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integrate  thermoelectric  phenomena  into  the  thermoelectric  mod¬ 
ule  as  well  as  heat  transfers  and  pressure  drops  in  the  heat 
exchangers.  First,  a  description  of  the  analytical  model  used  in  this 
paper  is  presented.  Based  on  this  model,  the  optimization  of  the 
heat  sink  design  (Dch,Nch)  and  the  operating  conditions 
4>h)  is  discussed  in  the  next  section. 

3.1.  Thermoelectric  phenomena 


In  1834,  Jean  Peltier,  a  French  physicist,  discovered  that  an  elec¬ 
trical  current  flowing  through  the  junctions  between  dissimilar 
materials  causes  a  temperature  difference.  Finally,  absorption  of 
thermal  power  </>c  at  the  cold  junction  and  generation  of  thermal 
power  </>h  at  the  hot  junction  are  observed.  These  two  thermal  pow¬ 
ers  are  defined  as  [15]: 

<t>h  =  a/r{  +  Ifiefecf2  -  K(TJh  -  T{)  (1) 


4>c  =  ccIT{  -  ^RdecI2  -  K(T[  -  T{)  (2) 

The  electrical  power  to  supply  at  the  thermoelectric  module  is 
the  difference  between  the  hot  and  cold  thermal  powers. 

4* elec  =  fill  ~  fie  =  ~  Vc)  +  RelecI  (3) 

The  performance  of  thermoelectric  modules  is  characterized  by 
the  COP,  which  is  defined  as  the  ratio  between  the  useful  thermal 
power  released  and  the  electrical  power  supplied.  In  heating  mode, 
the  useful  thermal  power  is  the  hot  thermal  power  and  the  source 
thermal  power  is  the  cold  thermal  power.  It  is  the  opposite  in  cool¬ 
ing  mode.  Finally,  depending  on  the  mode  considered  two  COPs 
can  be  written: 


COP™  =  -p-  =  A.  + 1  =  coptem  + 1  (4) 

4>elec  4)elec 

As  shown  by  Eqs.  (l)-(4),  the  performance  of  a  thermoelectric 
device  is  a  direct  function  of  the  material’s  properties  (a,  Reiec 
and  I<),  the  electrical  intensity  supplied  I  and  the  junction  temper¬ 
atures  TJh  and  Tjc.  The  material’s  properties  used  in  this  study  are 
provided  by  Ferrotec  [16]  and  correspond  to  bismuth  telluride: 

a  =  0.1 01  VIC1  (5) 

Reiec  =  1.7352  0  (6) 


I<  =  2.6637  W  K”1  (7) 

Fig.  3  plots  the  variation  of  COP  versus  the  hot  thermal  power 
demand,  for  different  temperature  differences  (with  TJC  =  285  K). 
It  clearly  shows  that  for  a  given  temperature  difference,  an  optimal 


Hot  thermal  power  (W) 

Fig.  3.  COP  as  a  function  of  the  hot  thermal  power  for  different  junction 
temperature  differences. 


hot  thermal  power  <f°hpt  exists,  leading  to  the  maximal  COP.  As  a 
consequence,  an  optimal  electrical  intensity  Iopt  exists,  leading  to 
the  maximal  COP.  The  greater  the  difference  in  junction  tempera¬ 
ture  is,  the  lower  the  maximal  COP  and  the  higher  the  correspond¬ 
ing  optimal  hot  thermal  power. 


3.2.  Heat  transfer  in  the  heat  exchangers 
3.2.1.  Heat  exchanger  design 

The  heat  exchangers  linked  to  the  hot  and  cold  junctions  are 
identical  to  ensure  the  reversibility  of  the  device  (heating  or  cool¬ 
ing  mode).  They  are  made  up  of  parallel  circular  mini-channels 
with  a  global  “Z”  configuration.  In  1996,  Kikas  [17]  demonstrated 
that  this  configuration  provides  a  good  distribution  of  the  thermal 
fluxes  on  the  entire  heat  exchanger  area.  In  fact,  a  good  distribution 
of  the  thermal  fluxes  provides  the  same  temperature  conditions  for 
all  the  semi-conductor  junctions.  Finally,  the  optimal  intensity 
supplied  to  the  TEM  is  optimal  for  all  the  semi-conductor 
junctions. 

The  calorific  fluid  used  is  water  and  the  thermal  conductivity  of 
the  heat  exchangers’  material  is  k  =  237  W  m_1  K_1  (aluminum). 
The  surface  area  (11)  occupied  by  the  mini-channels  corresponds 
to  the  thermoelectric  module  area  considered  (Fig.  4). 

Thus,  the  geometrical  constraints  enforce: 

D*=^l  (8) 

AT* 

The  height  H  corresponding  to  the  thickness  of  the  heat  exchan¬ 
ger  bulk  is  set  to  3  mm  for  this  study. 

One  of  the  weaknesses  of  this  kind  of  design  stems  from  the 
flow  distribution.  In  fact,  many  studies  have  shown  that  in  this 
configuration,  a  nonuniform  distribution  of  the  different  mass  flow 
in  the  mini-channel  is  observed  [18-21].  These  studies  illustrate 
that  the  greater  the  diameter  of  the  distributor  and  the  collector, 
the  better  the  flow  distribution  is.  We  can  also  note  that  the  distri¬ 
bution  depends  to  a  large  extent  on  the  number  of  mini-channels 
and  the  Reynolds  number. 

A  fluidic  model  was  built  and  integrated  into  the  global  model 
introduced  in  this  study  in  order  to  integrate  and  analyze  the  ef¬ 
fects  of  maldistribution  in  the  mini-channels  on  the  device’s  global 
COP.  Due  to  the  complexity  of  the  problem,  this  aspect  will  not  be 
presented  in  detail  in  this  paper.  However,  a  complementary  study 
has  shown  that  the  collector  and  distributor  diameters  must  be  set 
at  15  mm  to  obtain  uniform  distribution  of  the  fluid  in  all  the  mini¬ 
channels.  Consequently,  it  can  be  written: 


</c  = 


m 


h/c 


Nc 


(9) 


To  keep  the  same  temperature  difference  from  the  inlet  to  the 
outlet  of  both  heat  exchangers,  the  mass  flows  for  each  side  differ 
due  to  the  thermal  power  difference  between  the  hot  and  cold 
sides  (Eqs.  (1)  and  (2)). 

The  boundary  conditions  of  the  problem  are  presented  in  Fig.  5. 
The  thermal  fluxes  are  transferred  on  both  sides  of  the  central  heat 
exchanger.  Finally,  for  symmetry  reasons  the  study  domain  is  re¬ 
duced  to  one  mini-channel  on  its  half  height  for  the  hot  side. 


3.2.2.  Thermal  resistances 

The  total  thermal  resistances  R™c  from  the  junction  to  the  fluid 
can  be  divided  into  three  distinct  thermal  resistances  related  to  the 
thermal  conduction  through  the  bulk  of  the  heat  exchanger,  the 
convective  heat  transfer  in  the  mini-channels  and  the  heating  of 
the  fluid.  Finally,  the  total  thermal  resistances  for  the  hot  and  cold 
sides  are: 
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Fig.  4.  Heat  exchanger  design. 
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Fig.  5.  Boundary  conditions  of  the  problem. 


yj  yin 
)tot  _  J(1  ~  [h 

*  4>h 

yj  y  cer  ycer  y  b  y  b  yint  yint  y/  y f  yin 
=  1h~1h  +  Jb  -Jb  +  Jb-Jb  +  Jb  -lh  +  lh~lh 

$ h  4*h  4*1]  4*1]  4>h 


(10) 


R“f  =  Rcer  +  Rctc  +  Rcond  +  Rchonv  +  R{  (1 1 ) 

yOUt  yj 

^tot  _  J  c  J  c 
c  4>c 

•ycer  yj  yb  ycer  yint  yb  y/  yint  yOUt  y/ 

=  c  ~  c+  c~  c  +  c  r  c+  c~  c  +  c  r  c  (12) 

<Pc  (Pc  <Pc  <Pc  <Pc 

R™  =  Rcer  +  Kctc  +  Kcond  +  Rchonv  +  ^  (13) 

The  thermal  resistance  due  to  the  conduction  through  the  bulk 
of  the  heat  exchanger  is  estimated  as  follows: 

Rcond= ^  (14) 

The  thermal  constriction  resistances  due  to  the  circular  geometry 
of  the  channel  are  not  considered  in  this  study.  Consequently,  ther¬ 
mal  transfer  through  the  bulk  of  the  heat  sink  is  assumed  to  be  one¬ 
dimensional.  This  assumption  is  acceptable  because  the  constriction 
resistance  must  be  negligible  in  the  conduction  resistance  with  re¬ 
gard  to  heat  exchanger  dimensions  (H  =  6  mm  and  Dch  =  0.6  mm). 
Moreover,  we  show  that  the  convective  resistance  is  the  main  part 
of  the  total  thermal  resistance  (see  Fig.  10).  The  convective  thermal 
resistance  between  the  wall  of  the  mini-channel  and  the  fluid  flow  is 


different  for  the  heat  exchanger  through  which  the  cold  fluid  flows 
and  the  heat  exchanger  through  which  the  hot  fluid  flows.  In  fact, 
the  hot  heat  exchanger  at  the  center  of  the  device  has  to  dissipate 
two  thermal  fluxes  coming  from  its  two  opposite  sides.  Finally,  we 
can  consider  that  the  exchange  area  between  the  solid  and  the  fluid 
for  the  central  heat  exchanger  is  half  of  the  exchange  area  of  the  two 
heat  exchangers  at  the  extremity  of  the  device: 


nconv 


2 

Mint 


(15) 


j^conv  _ 


hcAir 


with  hh/c,  the  convection  coefficient  equals: 


hh/c  = 


Nuft/C7 

Dch 


(16) 


(17) 


and  the  Nusselt  number  (for  0.7  <  Pr  <  7  and  Re  <  2300)  equals  [22]: 

rB^V/3 


Nu/,/c  =  4.364  - 


0.086  (Reh/cPr^) 
l+0.1Pr(Reh/c^' 


(18) 


This  correlation  assumes  that  the  thermal  flux  is  circumferen¬ 
tially  uniform  around  the  mini-channel.  This  assumption  is  not 
fully  verified,  even  for  the  central  heat  sink.  A  complementary 
study  has  to  be  carried  out  to  evaluate  the  heat  transfer  coefficient 
precisely. 

However,  this  correlation  is  used  to  estimate  an  effective  aver¬ 
age  Nusselt  in  the  channel  in  this  paper.  The  optimal  results  may 
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be  slightly  different  estimating  the  convective  heat  transfer  coeffi¬ 
cient  differently,  but  the  approach  proposed  in  this  paper  remains 
valid. 

The  thermal  resistance  linked  to  the  heating  of  the  fluid  from 
the  inlet  T™/c  to  the  mean  fluid  temperature  in  the  channel  Tfh/C  is: 

^c  =  rh™Cp  (19) 


As  Tfh/C  is  the  mean  temperature  of  the  fluid  between  the  inlet 
and  the  outlet  of  the  mini-channel,  a  factor  of  2  appears. 

The  analytical  model  also  includes  the  thermal  resistances  com¬ 
ing  from  the  ceramics  of  the  TEM  and  the  contact  between  the  TEM 
and  the  bulk  of  the  heat  exchangers  from  [16,23]. 

Rcer  =  0.014 1<  vr1 

Rctc  =  0.003  KVT1 


3.2.3.  Pressure  drop 

Correlations  used  for  the  calculation  of  the  pressure  drops  are  a 
function  of  the  development  length,  defined  as  [24]: 


Ll”c  =  0.05R  eh/cDd 


(20) 


For  the  inlet  region,  we  have  (for  Re  <  2300)  [25]: 


.  13.74  «„)'«  +  + 


.+  \l/2 


0-5K/C 

with 


1  +0.00021  (x, 


h/cJ 


XZ/r  = 


l/Dc 


'h'c  Re,1/C  ^/c>'  and  xh/c  Re)i/c 

For  the  fully  developed  region,  we  have  [26,27]: 


APh/c  =J 


f  ((i-tf/M/cPf) 

Jh/C[  2  Dch  J 


(21) 


ifLZc<l  (22) 


(23) 


4.  Thermoelectric  heat  pump  optimization 


The  two  junction  temperatures  calculated  with  the  thermoelec¬ 
tric  model  are  linked  to  the  useful  and  source  fluid  temperatures 
(' Tfh  and  T{)  via  the  thermal  resistances  of  the  heat  exchangers. 


r(  -  t[ 

h  h 

n  tot 

(26) 

Tl  -  T{ 

c  c 

Rf 

(27) 

Fluid  temperatures  are  the  known  inputs  in  a  thermoelectric 
heat  pump  (rather  than  the  junction  temperatures,  which  are  most 
often  used  in  the  literature).  Eqs.  (26)  and  (27)  explicitly  show  that 
the  thermal  resistances  of  the  heat  exchangers  highly  influence  the 
junction  temperatures  and  therefore  the  device  performance.  In¬ 
deed,  Fig.  6  introduces  the  temperature  profiles  along  the  thermo¬ 
electric  heat  pump.  It  is  shown  that  the  lower  the  thermal 
resistance,  the  lower  the  temperature  difference  between  the  junc¬ 
tions  is,  resulting  in  increasing  the  COP  (see  Fig.  3). 

To  sum  up,  the  COP  of  a  TEM  integrated  into  a  thermoelectric 
heat  pump  is  a  direct  function  of  the  electrical  intensity  (see 
Fig.  1 )  supplied  and  of  the  heat  exchanger  thermal  resistances  link¬ 
ing  the  junction  temperature  to  the  fluid  temperature.  The  goal  of 
this  study  is  to  optimize  the  device’s  global  COP,  including  the 
pressure  drop  in  the  heat  exchangers.  The  device’s  COP  is  defined 
as: 


COP 


device  — 


0h 

c fielec  T"  0m 


(28) 


Finally,  to  optimize  the  device  we  need  to  optimize  both  the 
heat  exchanger  design  (Nch,  Dch )  and  the  operating  conditions 

( I°Pt ,  </c). 

First,  optimization  for  a  given  thermal  power  demand  is  han¬ 
dled  with  two  different  objective  functions  and  then  complete 
optimization  including  the  hot  thermal  power  as  a  complementary 
optimization  variable  of  the  device  is  presented.  The  direct  search 
method  [29]  is  used  with  the  Engineering  Equation  Solver  (EES) 
software. 


4A.  Optimization  for  a  fixed  thermal  power  demand 


with 


fh/c 


64 


(24) 


The  singular  pressure  drops  due  to  the  bifurcation  of  the  fluid 
from  the  distributor  to  the  channels  as  well  as  the  confluence  of 
the  fluid  from  the  channels  to  the  collector  are  not  treated  in  this 
study.  A  complementary  fluidic  model  including  ah  the  flow  distri¬ 
bution  phenomena  helps  estimate  the  error  caused  by  this 
assumption.  The  error  reported  is  about  20%,  so  that  singular  pres¬ 
sure  drops  are  ignored  to  reduce  the  complexity  of  the  model.  The 
authors  are  aware  that  the  error  induced  by  this  assumption  may 
modify  the  optimal  design.  However,  the  increase  in  total  pressure 
drop  influences  the  optimal  geometry  of  the  heat  exchanger  very 
slightly,  as  reported  in  [28]. 

Finally,  the  total  mechanical  consumption  of  the  two  fluid  circu¬ 
lating  auxiliaries  (assuming  a  pump  efficiency  of  1 )  corresponds  to: 


^  mfAP,  +  2mfAPc 

0m  =  — - - -  (25) 

As  described  in  Fig.  1,  there  is  one  heat  exchanger  at  the  center 
of  the  device  through  which  the  hot  fluid  flows  and  the  two  heat 
exchangers  at  the  extremity  through  which  the  cold  fluid  flows. 


First  let  us  consider  a  case  with  a  fixed  thermal  power  demand. 
The  optimization  problem  can  be  handled  by  maximizing  the  de¬ 
vice’s  COP.  The  optimization  variables  are  the  mass  flows  for  the 
cold  and  hot  sides,  the  channel  diameter  and  the  number  of 
channels. 

max  COPdeM-ce (X)  e  U  (29) 

XeR4 


P^CCCCS 


Fig.  6.  Influence  of  the  thermal  resistances  on  the  junction  temperatures. 
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X  =  (Dc\  Nch,  mffc) 

(30) 

Subjected  to: 

E  >  0.0002  m 

(31) 

Dch  >  0.0002  m 

(32) 

< 

II 

"P? 

< 

(33) 

Parameters  E  and  Dch  are  constraints  due  to  the  limit  of  the 
manufacturing  process.  The  optimization  results  show  that  what¬ 
ever  the  useful  thermal  power,  the  optimal  value  of  E  corresponds 
to  the  lower  bound.  A  constant  temperature  difference  between 
the  hot  and  cold  sides  is  assumed,  as  has  been  explained  in  Section 
3.2.1.  Fig.  7  shows  that  the  optimal  number  of  channels  decreases 
with  the  increase  of  the  hot  thermal  power  demand.  As  a  conse¬ 
quence,  and  due  to  the  fact  that  E  is  constant  and  equal  to  the  low¬ 
er  bound  value,  the  mini-channel  diameter  and  finally  the  total 
exchange  area  between  solid  and  fluid  increase  with  the  hot  ther¬ 
mal  power  demand  (see  Eq.  (8)).  The  optimal  mass  flow  in  each 
heat  exchanger  increases  with  the  increase  of  the  hot  thermal 
power  demand  in  order  to  keep  the  heat  transfer  coefficient  high 
(Fig.  8).  Finally,  an  increase  of  the  optimal  Reynolds  number 
(Fig.  8)  and  thus  the  pressure  drop  and  the  heat  transfer  coefficient 
is  observed  when  the  hot  thermal  power  increases  (Fig.  9). 

Fig.  10  shows  the  relative  proportion  of  the  convection  resis¬ 
tance  to  the  conduction  resistance  as  a  function  of  the  thermal 
power  demand.  The  higher  the  thermal  power  demand,  the  lower 
this  ratio  is.  Indeed,  the  thermal  resistance  due  to  conduction  is 
constant  but  the  higher  the  thermal  power,  the  higher  the  convec¬ 
tive  heat  transfer  coefficient  is.  The  ratio  is  within  the  range  of  5-8. 
For  this  reason,  the  controlling  thermal  resistance  is  the  convec¬ 
tion,  as  stated  in  Section  3.2.2,  and  the  constriction  resistance  in 
the  bulk  can  be  ignored. 

Fig.  1 1  clearly  illustrates  the  weak  influence  of  the  mechanical 
consumption  of  the  fluid-circulating  auxiliaries  compared  to  the 
electrical  consumption  of  the  thermoelectric  modules.  It  helps  to 
validate  the  assumption  related  to  ignoring  the  singular  pressure 
drops,  stated  in  Section  3.2.3. 

Fig.  12  illustrates  the  optimized  COP  versus  hot  thermal  power 
demand.  The  optimal  COP  shows  a  maximum  corresponding  to  an 
optimal  hot  thermal  power  demand  of  about  55  W  (for  the  case 
considered).  The  changes  in  COP  versus  hot  thermal  power 
demand  corresponds  to  the  classical  behavior  of  a  thermoelectric 
element  versus  hot  thermal  power,  as  shown  in  Fig.  3. 

One  of  the  most  widely  used  criteria  for  the  optimization  of  the 
heat  exchanger  [30,31]  and  for  the  thermoelectric  devices  [32,33] 


Fig.  8.  Optimal  mass  flow  for  the  cold  and  hot  sides  as  a  function  of  the  hot  thermal 
power  demand. 


Fig.  9.  Heat  transfer  coefficient  and  Reynolds  number  as  a  function  of  hot  thermal 
power  demand. 
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Fig.  10.  Ratio  between  the  thermal  resistance  due  to  convection  and  conduction  as 
a  function  of  hot  thermal  power  demand. 


Fig.  7.  Optimal  number  of  channels  and  optimal  channel  diameter  as  a  function  of  Fig.  11.  Ratio  between  electrical  consumption  and  mechanical  consumption  as  a 
the  hot  thermal  power  demand.  function  of  hot  thermal  power  demand. 
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Fig.  12.  Optimal  COP  of  the  device  as  a  function  of  the  hot  thermal  power  demand. 

Fig.  14.  Optimal  COP  of  the  device  as  a  function  of  the  hot  thermal  power  demand 
and  for  two  optimization  criteria. 


has  long  been  entropy  generation  minimization  (EGM).  For  this 
study,  an  adiabatic  device  is  considered,  so  the  two  principal  entro¬ 
py  generation  sources  result  from  the  pressure  drop  in  the  heat 
exchangers  and  the  thermal  transfer.  The  three  subsystems  gener¬ 
ating  entropy  are  described  in  Fig.  13. 

In  the  thermoelectric  element,  the  entropy  generation  due  to 
thermal  transfer  is: 


<-TEM  _  _  (frc  , 

Sen  ~  jj  +  p 


(34) 


For  the  cold  and  hot  sides  of  the  heat  exchangers,  the  entropy 
generation  is,  respectively  [34]: 


1  - rOUt 

- f  +  mtcotCpLn-^-  + 
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m'otAPc 
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(36) 


Finally  the  total  entropy  generation  of  the  device  is  the  sum  of 
the  entropy  generation  of  each  subsystem,  which  reduces  to: 


mt°tCpLn^-  + 


mg*APh 

pTh 


■  mc°cCpLn  + 


m'otAPc 


(37) 


of  the  results  between  these  two  optimization  criteria  leads  to 
concluding  that: 

max  COPdewce(X)  e  R  «  COP(minS™,(X)  e  R)  (40) 

XeK4  XeR4  & 

Compared  to  the  optimal  configuration  corresponding  to  the 
maximization  of  the  device’s  COP  and  the  EGM,  a  slight  difference 
is  observed,  mainly  for  the  optimal  number  of  channels  (an 
approximately  8%  difference)  (Figs.  15  and  16).  In  fact,  Fig.  17 
shows  that  for  a  hot  thermal  power  set  to  55  W  the  range  of  chan¬ 
nel  numbers  leading  to  the  optimal  COP  is  quite  wide.  This  low 
sensitivity  of  the  maximal  COP  to  the  number  of  channels  explains 
that  it  is  possible  to  have  a  variation  of  the  optimal  configuration 
between  the  two  optimization  methods  without  substantially 
varying  the  device’s  resulting  optimal  COP. 

The  optimization  for  given  hot  thermal  power  shows  that  there 
exists  an  optimal  configuration  of  the  heat  exchangers  and  optimal 
mass  flow  for  the  hot  and  cold  sides.  The  results  also  show  that 
there  is  an  optimal  hot  thermal  power  leading  to  the  maximal  opti¬ 
mal  COP  of  the  device.  This  last  observation  turns  us  to  the  follow¬ 
ing  complete  optimization,  including  the  hot  thermal  power  as  a 
complementary  variable. 

4.2.  Complete  optimization 


The  optimization  problem  can  thus  be  written: 
minS“(X)  e  R  (38) 

XeR4 

X=(Dch,  Nch,  mc°c,  m‘hot)  (39) 


In  this  part,  the  hot  thermal  power  demand  is  considered  as  an 
optimization  variable  in  order  to  validate  the  optimal  configuration 
obtained  by  partial  optimization  (29).  Finally  the  optimization 
problem  can  be  written: 


The  constraints  are  the  same  as  for  the  maximization  of  the 
device’s  COP  (Eqs.  (31)-(33)). 

Fig.  14  compares  the  device’s  COP  obtained  by  maximizing 
the  device’s  COP  and  by  using  the  EGM  criterion.  The  similarity 


Fig.  13.  Decomposition  of  entropy  generation  sources. 


-a-  Optimal  number  of  channels  with  COP  maximization 
Optimal  number  of  channels  with  EGM 
-a-  Optimal  channel  diameter  with  COP  maximisation 
Optimal  channel  diameter  with  EGM 


Fig.  15.  Comparison  of  optimal  geometries  obtained  by  maximizing  the  device’s 
COP  and  by  EGM. 
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-o-  Optimal  mass  flow  for  the  hot  side  with  COP  maximization 
Optimal  mass  flow  for  the  hot  side  with  EGM 
-a-  Optimal  mass  flow  for  the  cold  side  with  COP  maximization 


20  40  60  80  100 

Hot  thermal  power  demand  (W) 

Fig.  16.  Comparison  of  optimal  mass  flows  obtained  by  maximizing  the  device’s 
COP  and  by  EGM. 


Number  of  channels 

Fig.  17.  Optimal  COP  of  the  device  as  a  function  of  the  number  of  channels  for  a 
thermal  power  demand  of  55  W. 

max  COP^ice  (X)  e  U  (41) 

XeR5 

X  =  (Dc\  Nch,  mf,  m‘hot,  <t>h)  (42) 

The  constraints  are  the  same  as  for  the  previous  cases  ((33)). 
For  this  optimization,  the  following  results  are  obtained: 

Doptch  =  6.210-4  m;  Nopt  ch  =  54;  =  0.02629  Kg  s"1; 

m'ot;opt  =  0.01562  Kg  s"1 ;  fhpt  =  54.31  W  giving  : 

C0p0device  =  2  A4:  RT  =  0  04  K  W_1 ;  RT  =  0  02  K  W1 ; 

A P„  =  8050  Pa;  A Pc  =  6170  Pa 

The  optimal  thermal  power,  the  mass  flows,  the  diameter  and 
the  number  of  mini-channels  obtained  with  the  complete  optimi¬ 
zation  is  the  same  as  those  obtained  with  the  partial  optimization. 

Please  note  that,  with  complete  optimization,  including  the 
thermal  power  as  an  optimization  variable,  the  EGM  criterion  can¬ 
not  be  used.  In  fact,  the  minimum  entropy  generation  of  a  thermal 
device  is  reached  when  the  thermal  power  is  minimal.  Thus, 
according  to  the  EGM  criterion,  the  optimal  operating  condition 
is  reached  when  the  thermal  power  is  null. 

5.  Conclusion 

Thermoelectric  heat  pumps  are  an  innovative  technology  for 
heating  or  cooling  and  show  a  high  potential  for  low-consumption 
buildings.  High  thermal  power  can  easily  be  reached  by  parallel 
association  of  the  optimized  unit  presented  in  this  paper.  One  of 


the  most  critical  aspects  of  these  devices  is  the  many  thermal  resis¬ 
tances  from  the  junctions  of  semi-conductors  to  the  working  fluid. 

In  this  paper,  an  original  configuration  of  a  thermoelectric  heat 
pump  is  introduced.  The  originality  of  this  configuration  comes 
from  its  compactness  due  to  a  thermal  transfer  on  both  sides  of 
the  central  heat  exchanger  and  to  its  modularity.  In  fact,  this  con¬ 
figuration  provides  two  operating  modes:  TEM  association  in  par¬ 
allel  or  TEM  association  in  cascade. 

This  paper  focuses  on  the  parallel  operating  mode.  An  analytical 
model  was  developed  in  order  to  completely  optimize  the  device. 
The  optimization  criteria  considered  is  the  global  COP  including 
electrical  consumption  of  the  circulating  auxiliaries. 

The  first  part  of  the  optimization  introduced  in  this  study  high¬ 
lights  that  an  optimal  heat  exchanger  design  exists  allowing  one  to 
maximize  the  device’s  global  COP.  In  fact,  when  the  hot  thermal 
power  demand  increases,  the  optimal  exchange  area  between  the 
solid  and  the  fluid  increases  as  well.  This  increase  is  due  to  a  de¬ 
crease  of  the  number  of  channels  and  an  increase  of  their  diameter 
due  to  a  constant  space  between  them.  Concerning  the  optimal 
mass  flow,  an  increase  of  its  value  is  observed  with  the  increase 
of  the  hot  thermal  power  demand.  Finally,  when  the  hot  thermal 
power  increases,  the  optimal  design  (number  of  channels)  and 
operating  conditions  associated  (mass  flow)  increase  the  Reynolds 
number  and  thus  the  pressure  drop  and  the  heat  transfer 
coefficient. 

These  results  also  show  that  for  given  fluid  temperature  condi¬ 
tions,  there  exists  one  optimal  thermal  power  leading  to  the  max¬ 
imal  COP. 

The  second  result  introduced  in  this  paper  is  that  it  is  possible 
to  optimize  the  device  by  EGM.  In  fact,  by  minimizing  this  crite¬ 
rion,  the  same  COPs  are  obtained  by  maximizing  the  device’s  global 
COP  and  by  EGM.  There  is  little  variation  of  results  concerning  the 
corresponding  optimal  design  due  to  the  low  sensitivity  of  the 
optimal  COP  to  the  number  of  channels  within  a  large  range. 

Finally,  it  is  possible  to  completely  optimize  the  device,  includ¬ 
ing  the  hot  thermal  power  as  an  optimization  variable.  Conse¬ 
quently,  we  obtained  an  optimized  device  for  given  temperature 
conditions  with  the  possibility  of  extending  the  thermal  power 
generated  with  an  association  in  parallel  or  in  cascade  without 
degrading  the  optimal  COP. 

In  a  future  perspective,  a  flow  distribution  model  must  be  added 
to  take  into  account  the  pressure  drop  caused  by  the  collection/dis¬ 
tribution  of  water  in  the  mini-channels.  The  optimized  results  may 
change  with  regard  to  the  total  pressure  drop,  including  the  collec¬ 
tor/distributor,  but  the  change  may  be  slight  as  the  mechanical 
consumption  of  the  auxiliaries  is  two  orders  of  magnitude  below 
the  electrical  consumption  of  the  thermoelectric  modules. 

Similarly,  the  estimation  of  the  heat  transfer  coefficient  must  be 
improved.  Optimized  results  may  differ  slightly. 
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